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We perform a linear stability analysis of dynamical Chern-Simons modified gravity in the geo- 
metric optics approximation and find that it is linearly stable on the backgrounds considered. Our 
analysis also reveals that gravitational waves in the modified theory travel at the speed of light in 
Minkowski spacetime. However, on a Schwarzschild background the characteristic speed of propaga- 
tion along a given direction splits into two modes, one subluminal and one superluminal. The width 
of the splitting depends on the azimuthal components of the propagation vector, is linearly propor- 
tional to the mass of the black hole, and decreases with the third inverse power of the distance from 
the black hole. Radial propagation is unaffected, implying that as probed by gravitational waves the 
location of the event horizon of the spacetime is unaltered. The analysis further reveals that when a 
high frequency, pure gravitational wave is scattered from a black hole, a scalar wave of comparable 
amplitude is excited, and vice-versa. 



Introduction. The unknown nature of dark matter and 
dark energy, and the possibility that one or both of these 
phenomena could be a misinterpretation of observations 
if general relativity is not the correct theory of spacetime 
on all scales, has renewed interest in alternative theo- 
ries of gravity. One particular alternative theory that 
has attracted considerable attention recently is dynam- 
ical Chern-Simons (CS) modified gravity [l|. This the- 
ory is interesting because it passes all Solar System and 
binary pulsar constraints, while allowing for large devi- 
ations in highly curved spacetimes, such as for exam- 
ple around rotating black holes [2, Isj. Moreover, terms 
akin to the CS modification in the action generically arise 
in string theory [J, |5| , effective field theory [6[ and also 
possibly in loop quantum gravity [7|, |8|. In this sense, 
CS modified gravity can be considered theoretically well- 
motivated. The theory is also an excellent proxy for the 
study of gravitational parity violation, as it includes the 
only non- vanishing curvature-squared term that violates 
this symmetry 1]. 

Chern-Simons modified gravity extends the Einstein- 
Hilbert action through the addition of the product of a 
scalar field with the Pontryagin invariant, which is the 
contraction of the Riemann tensor and its dual. In the 
dynamical theory studied here, the action also includes a 
kinetic term for the scalar field. The Pontryagin invari- 
ant is a curvature squared correction that breaks parity 
and time-reversal invariance, due to the Levi-Civita ten- 
sor in the definition of the Riemann dual. Such a term 
leads to modifications to the field equations that alter 
General Relativity (GR) solutions that break parity. Ex- 
amples of such GR solutions are the Kerr spacetime or a 
binary compact object spacetime, where either the spin 
or the total angular momentum vectors select a preferred 
direction that couples to the CS correction. Thus, for ex- 



ample, though Schwarzschild is a solution to CS modified 
gravity, Kerr is not. 

The analytical complexity of the CS modified field 
equations makes it almost impossible to find closed- 
form solutions to astrophysically interesting spacetimes. 
An alternative is to evolve such spacetimes numerically. 
However, the modified field equations generically include 
third time derivatives in the form of contractions of the 
Levi-Civita tensor with the covariant derivative of the 
Ricci scalar, which could render the theory ill-posed and 
unstable. In fact, Ostrogradski showed many years ago 
that theories defined by Lagrangians with more than one 
time derivative, which cannot be eliminated by partial 
integration, are linearly unstable |9]]. If one considers CS 
modified gravity to be the leading order correction to 
GR of some putative well behaved gravitational theory, 
then it might be reasonable to expect that the solutions 
of the stable sector of the effective theory are good ap- 
proximations to the physically-relevant solutions of the 
full theory. In principle, the stable sector can be found 
by order reduction, i.e., performing a perturbative expan- 
sion in the coupling parameter of the effective theory jlO| . 
A numerical evolution of CS modified gravity will there- 
fore need to implement a similar algorithm, rather than 
attempt a "brute-force" solution. 

Whether dynamical CS gravity, even when considered 
as an effective theory, is mathematically well-posed is un- 
known at present. One necessary (though not sufficient) 
condition for well-posedness is linear stability of the the- 
ory about backgrounds that are stable in GR. Clearly, if 
the background is fiat, then CS gravity is linearly stable 
as higher derivatives can only arise at higher than lin- 
ear order, since the Riemann tensor vanishes. Thus, a 
meaningful test of linear stability requires study of per- 
turbations about a solution to CS modified gravity that 



has non-trivial curvature. 

In this paper we perform a Unear stabihty analysis 
in the geometric optics (high-frequency) approximation. 
We derive general relationships valid for an arbitrary 
background geometry, then evaluate these in a few spe- 
cific spacetimes, including Minkowski and Schwarzschild. 
We show explicitly that dynamical CS gravity is linearly 
stable on these backgrounds, though for Schwarzschild 
this is only true when the sign of the coupling constant 
is chosen such that the CS scalar field possesses positive 
energy^. As a byproduct of this analysis, we derive the 
speed of coupled gravitational wave/CS-scalar modes, 
and show that, while in Minkowski spacetime it is exactly 
the speed of light, in a Schwarzschild spacetime it differs 
from this speed, allowing for both subluminal and super- 
luminal modes. One way to physically understand this 
is to consider the CS-modifications to the Einstein field 
equations as an additional effective stress-energy tensor 
for the gravitational field; this effective matter violates 
the usual energy conditions, allowing for superluminal 
modes. 

In this paper we follow mostly the conventions of [12| . 
We use abstract index notation, the metric signature 
(— ,-|-,-f,-|-), geometric units where G = c = 1, and 
overhead bars to denote quantities associated with the 
background spacetime. 

Basics of CS gravity. The CS modified action leads to 
the following field equations (we refer the reader to [l] 
for the action and other details): 
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where Gab is the Einstein tensor, D is the D'Alembertian 
operator, k = (IOttG)"^, and the C-tensor is 
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with the dual Riemann tensor defined as 

1 
2' 

and e^'^'^/ the 4-dimensional Levi-Civita tensor. In this 
study we do not consider any potential for the CS scalar 
(see e.g., jl[). The Pontryagin invariant is defined as 

*RR = *R"'b'^ R acd — -Z(^ '^ R°'befR acd , (5) 



while the total stress-energy tensor is 

rp rpmat 1^ rp^ 

J^ab — -l^ab + -<- ab^ 



(6) 



^ Incidentally, this physically-reasonable choice is opposite to 
the sign used in an earlier study of quasi-normal modes on a 
Schwarzschild background [llj . which explains why [lll | found 
unstable modes 



where T^'^* is the matter stress-energy, which we set to 
zero here [1[ , and T^j, is the stress-energy of the CS field 
I?, given by 
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The quantities (a, /3) are coupling constants whose di- 
mensions depend on the dimensions one chooses for the 
scalar field [l|, |2| . We shall find it convenient to define 
the parameter ^ = a^ /{(3k), as CS corrections will be 
proportional to it. 

If we move the C-tensor to the right-hand-side of the 
modified field equations, and consider this an effective 
stress-energy tensor, then the energy conditions are not 
necessarily satisfied in an arbitrary background. One 
can see that the sign of x'^X^Gab depends not only on 
x'x'T^arS but also on x^x'T^b and x'^x'Cab, where x" 
is a timelike or null vector, depending on the energy 
condition considered. In certain non-trivial spacetimes, 
such as the spinning black hole solution recently found 
in [2, ISj or the perturbed Schwarzschild spacetime con- 
sidered here, the energy conditions are not satisfied. 

Linear Stability Analysis. We study the modified field 
equations in perturbation theory, decomposing the full, 
spacetime metric into gab = gab + e/i hab and the full CS 
scalar field into ■& = d + eh S-d. We work to first order 
in the perturbation book-keeping parameter e^, which 
labels the order of the perturbation. The background 
quantities {gab,^) sue solutions to the unperturbed field 
equations, while (hab, Si!}) are small perturbations away 
from this background. 

We consider background solutions where the metric is 
vacuum and where both the scalar field and the Pontrya- 
gin invariant vanish. Examples of such solutions are the 
Minkowski and Schwarzschild spacetimes. We seek solu- 
tions to the perturbed metric and CS field of the form 



5i9 = B(t,x^)e''^(*'^'°)/<^* 
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We further impose the geometric optics approximation, 
where we require that the phase (j) varies much faster than 
the amplitudes (Aab^B). This is enforced by requiring 
that the phase book-keeping parameter e^ ^ 1. With 
the above ansatz we are also restricting this analysis to 
steady-state solutions, hence the CS scalar and metric 
perturbation have the same phase. 

The modified field equations and equation of motion 
for the CS field can now be expanded bivariately in e^ ^ 
1 and £0^1. The dominant term in the expansion is of 
0{eh, e|). The perturbed equation of motion for the CS 
field becomes 
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while the perturbed CS field equations become 

habk'k, = -{Sd)e''^f^aCb)'efkgkd . (10) 



Here we have defined the four wavevector k"- = V°0, and 
Cabcd is the background Weyl tensor (which is equal to 
the background Rieniann tensor since the background is 
vacuum) . 

The perturbed field equations and the equation of mo- 
tion for the perturbed scalar field can be combined to 
obtained a modified dispersion relation, namely 

(Kq/C j = ~^£ £ d^abc ^f ijf^ f^ekgiih 

= 2i {*CfJ) {*Chdfg) k-k^k^k^ . (11) 

For notational convenience, we define the tensors Sab and 
Wab in analogy with the magnetic and electric parts of 
the Weyl tensor 
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Notice that Wac plays the role of the electric part of the 
Weyl tensor Eab and Sac that of the magnetic part Bab, 
except that here the Weyl tensor is contracted with the 
four- wavevector A:° , instead of a timelike unit vector [13| . 
In terms of these quantities, the modified dispersion re- 
lation can be written in the following simple form 



[kak'^y = 2^ iS'"') (Sab) 



(13) 



Some interesting conclusions follow from these results. 
First, the right-hand side of Eq. (fT3|) vanishes trivially 
for all Petrov type-0 (conformally flat) spacetimes, such 
as Minkowski, simply because the Weyl tensor vanishes. 
Therefore, GWs in dynamical CS gravity have the usual 
dispersion relation fcafc" = to linear order in these 
spacetimes and they travel at exactly the speed of light. 
Second, when ^ > 0, which is natural as this would im- 
part positive energy to the scalar fleld, then linear stabil- 
ity requires that S'"^ S^ > on the given background. 
Using the results of [l3|, it is straight foward to show 
that all Petrov type-N and type-Ill (pure gravitational 
wave) spacetimes have kak"' > to linear order. If not 
strictly 0, wave speeds will differ from that of light. For 
spacetimes that have less symmetries in the Weyl tensor, 
the right-hand side of Eq. (fT3|) does not vanish in general, 
and we are not aware of results regarding its positivity. 
Such spacetimes would need to be examined on a case- 
by-case basis, as we do next for Schwarzschild. 

The Speed of Gravitational Waves in Dynamical CS 
Theory on a Schwarzschild Background. The modified 
dispersion relation of Eq. (fTB]) can be solved given a spe- 
cific background that is a solution to CS modified gravity, 
and here we do so for the Schwarzschild metric. Param- 
eterizing the four wave- vector as fc" = [17, /ci, /c2, fcs], we 



find that in a Schwarzschild background, keeping only 
the correction that is lowest order in C, 
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(14) 



Here (i, r, 9, 0) are standard Schwarzschild coordinates, 
m is the black h( 
the GR result is 



m is the black hole mass, f — I — 2m/r, ( — S^/rn^ and 
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This solution to the CS modified dispersion relation 
is different from the amplitude birefringent one found 
in [ij, ll5| , because here we are considering the dynami- 
cal theory, instead of the non-dynamical one and we are 
treating Si) as a perturbation, with no background ^. 

This solution deserves further discussion. First, notice 
that when dealing with an ingoing or outgoing spherical 
wavefront, such that A;2,3 — 0, there is no CS correction, 
n — r^GR. This makes sense, as a spherical wave front 
respects the spherical symmetry of the background, not 
providing the CS correction with any vector to couple to. 
Conversely, the maximal deviation from GR is achieved 
when there is no radial part to the wave vector. Second, 
the correction vanishes asymptotically and at the event 
horizon r — 2m (/ — 0), as can be seen by substituting 
Eq. ([T5|) into Eq. (|14|) . Thus, gravitational waves still 
"see" the same event horizon of the underlying geometry. 
Third, the CS correction to ficR is a term of third post- 
Newtonian order, as it is proportional to [Gm/{c'^r)]^, 
momentarily reinstating physical units. Fourth, the cor- 
rection scales as C^^^: instead of ( as found for the CS 
modified Kerr solution in [2], implying the effect found 
in this paper is, in this sense, larger. 

Perhaps the most interesting feature of the solution 
found above is the appearance of sub- and superlumi- 
nal modes. The modified dispersion relation is a quartic 
equation for f], instead of a quadratic one as in GR. Thus, 
one now finds four independent propagating modes: two 
representing ingoing and outgoing waves, and for each 
of these, two more representing subluminal and super- 
luminal propagation. In fact, one can check that the 
characteristics are no longer null, but now satisfy 
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The most general solution is a superposition of all 
modes: 

"ab — ^ab ^ ^ ^ab ^ ^ ^ab ^ ^ ^ab ^ ' 

S^ = B°'+e"t'+ + B°'-e''l'- + S''+e"*'^+ + S^'^e"'"^- , 

(17) 

where the phase (j)± is subject to Eq. P^ with either 
the plus or minus sign. In the regime where k°'ka ^ 



the amplitude coefficients are not independent quantities, 
but instead related via Eqs. ((9|)- ((T0)) . which reduce to 



A 2a i? 






because SabS''^ = {kck'=)^ / {2^) by Eq 1^. Part of the 
reason for this tight coupling between the fields is that 
we are only looking at steady-state solutions. In other 
words, nothing prevents us from choosing initial condi- 
tions to the original system ([T]) that consists of an arbi- 
trary combination of CS scalar perturbation plus gravita- 
tional wave perturbation. However, due to the coupling 
in the evolution equations, gravitational wave perturba- 
tions will source CS scalar perturbations, and vice- versa, 
and if eventually a steady state is reached, the ratios of 
amplitudes will be given by Eq. (|18l) . This implies, for ex- 
ample, that the scattering of a pure gravitational wave off 
a Schwarzschild black hole in CS modified gravity would 
source an outgoing scalar field component. 

Summary. In this paper we have shown that dynami- 
cal CS gravity is linearly stable to high-frequency pertur- 
bations on a variety of physically relevant backgrounds. 
An interesting further consequence, valid on general non- 
trivial backgrounds, is that the speed of gravitational 
waves is not equal to that of light. One might worry that 
such superluminal propagation violates causality. Super- 
luminal propagation of a field, however, is not necessarily 
related to causality violation. In particular, the perturba- 
tion to the metric is small, and so the causal structure of 
the spacetime is essentially unaltered. Furthermore, one 
cannot use the gravitational field to communicate signals 
acausally — that would require arbitrarily fast speeds of 
propagation, which is not the case here (see for exam- 
ple [16| for a similar discussion on such issues arising in 
k- essence cosmological models). 

One possible effective interpretation is to think of the 
background geometry as definiing a preferred reference 
frame, or "aether", in which gravitational wave distur- 
bances propagate, breaking local Lorentz invariance. In 
this sense, one can think of CS gravity as promoting 
the background spacetime to a non-vacuum medium on 
which gravitational waves propagate with an anomalous 
dispersion relation, similar to the analogous situation 
of electromagnetic waves propagation in an anomalous 
medium. 

An interesting consequence of such anomalous GW dis- 
persion is the plausible temporal separation of GW and 
electromagnetic pulses that would otherwise be coinci- 
dent in GR. Our results suggest that as a GW scatters 
off a non-trivial background, two modes will generically 
be excited (a super-luminal and a sub-luminal one) that 
will separate as these waves traverse the non-trivial back- 
ground. For example, if the computed CS effect per- 
sists when galaxies act as the lens of cosmological GWs, 



then these waves might accumulate a delay. However, 
a back of the envelope estimate, together with current 
constraints on $,-^/^ < 10"* km [2|, suggest that this ef- 
fect might not currently be measurable. This is because 
C — ^/m^ ^ 1 when one saturates the current constraints 
on ^ and uses galactic masses on the order of 10^"^ Mq. 
As argued in |2|, |l7 1 , the best constraints on dynamical 



CS gravity arise for small compact objects, where the 
Riemann tensor, which scales as the inverse of the mass 
squared, is large. 
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